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The results of extensive microscopic Relativistic Mean Field (RMF) calculations for the 
nuclei appearing in the a - decay chains of recently discovered superheavy elements with 
109 < Z < 118 are presented and discussed. The calculated ground state properties like 
total binding energies, Q values, deformations, radii and densities closely agree with the 
corresponding experimental data, where available. The root mean square radii closely fol- 
low A 1 / 3 law (A being the mass number) with the constant r = 0.9639 ± 0.0005 fm. The 
double folding (tpp) approximation is used to calculate the interaction potential between 
the daughter and the a, using RMF densities along with the density dependent nucleon - 
nucleon interaction (M3Y). This in turn, is employed within the WKB approximation to 
estimate the half lives without any additional parameter for a - decay. The half lives are 
highly sensitive to the Q values used and qualitatively agree with the corresponding ex- 
perimental values. The use of experimental Q values in the WKB approximation improves 
the agreement with the experiment, indicating that the resulting interaction potential is 
reliable and can be used with confidence as the real part of the optical potential in other 
scattering and reaction processes. 

1. Introduction 

The synthesis of the superheavy elements (SHE) is a subject that has remained at 
the forefront of nuclear physics ever since the prediction of islands of stability in this 
mass region (Z ~ 114, N ~ 184) in the 1960's [d I2]- After about four decades of 
experimentation, it still does not seem to be possible to reach the N = 184 closed shell 
with stable neutron rich beams such as 48 Ca used at JINR/Dubna. In order to get closer, 
radioactive ion beams with neutron excess higher than 48 Ca will be required such as 
those that may become available with the DRIBs facility at Dubna [ 4] and elsewhere. 

We start with a short review of the current status of experimental work in Section II, 
followed by a brief description of the theoretical framework in Section III. The results 
and discussion are contained in section IV. The conclusions are presented in the last (V) 
section. 



2. Status of Experimental Work 

The elements above Fermium have been created in the laboratory by completely fus- 
ing heavy ions. The evaporation residues (EVR's) of compound nuclei (CN) upto about 
Z = 112 are produced primarily by successive alpha - decays to known nuclei enabling 
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their unambiguous identification by the method of a - a correlations [ |3] . However, at 
Seaborgium, the production cross - sections fall below one nanobarn with half - lives of 
less than a second. To access the new region above Z=106, new methods had to be devised 
[E1IZ1 such as 'cold-fusion' using Pb and Bi targets [EJ, invented by the Dubna Group and 
successfully used in the early discoveries of Bh, Hs and Mt. Together with increasingly 
sensitive detection methods it thus became possible to identify a single atom [EJ. Using 
such techniques, two isotopes of Hassium, 264 > 265 Hs were created at GSI in the 1980's [ 
FTTH ITT] . 267 Hs in the mid - 90's at Dubna [UJ and 277 Hs as a part of the Z = 114 chain 
[H3j also at Dubna. The GSI group had successfully produced the elements 270 Ds [Ej, 
269,27i DS; 272 111; [ [j^j [jjj [j^j [jj by 1998 and 277 m [[jjj by the year 2 002. Additionally, 

two decay chains have come from LBNL confirming the 271 Ds data from GSI using the 
same 208 Pb( 64 Ni,n) reaction [|20] but with gas-filled detectors to enhance EVR collection. 
In the last couple of years new results have become available from RIKEN for 271 110 [ 
I2TJI22], 272 1H [123 and for 278 113 [ 24 . With extensive activity underway at this facility, 
the current experimental program includes the on - going cold fusion synthesis of recently 
completed experiments for 277 112 and future plans for 283 114 Thus, cold fusion 

reactions usually involving the evaporation of one neutron, have successfully been used to 
synthesize elements up to Z=112. However the maximum production cross-section for this 
element (Z = 112) is achieved at an excitation energy of only ~ 10 MeV. If one were to ex- 
trapolate, to reach the heavier elements the excitation energy at the Coulomb barrier may 
well drop below zero (see for an exhaustive discussion). Clearly, as one approaches 
the predicted island of stability at N=184 (assuming the correctness of this postulate), 
a higher excitation energy is required. This is achievable through greater neutron excess 
in both the target and the projectiles as with 'hot fusion' reactions. The Dubna group 
have successfully used actinide targets with a 48 Ca beam to fuse Z>114, where the com- 
pound nucleus is formed following the evaporation of three or more neutrons. Elements 

286,287,288,289 114; 287,288 -Qg and 290,291,293 -Qg have been produced at higher excitation 

energies and different cross-sections [|23l2Ell2HlEn]- Preliminary results for 294 118 from a 
single event have also been cited [EH] with further experimentation in progress. As a final 
observation, even with hot fusion reactions, EVR's may not pass through a region greater 
than N=174-177 [E2]- Though this is closer by about 6-8 neutrons to the N=184 region, 
it may not be sufficient. It is hoped that with the dawn of the new era of radioactive ion 
beams the 'sea of instability' may be breached. 

From the above discussion, it is clear that above Z=112, especially in reactions involving 
neutron rich beams of the rare 48 Ca isotope, a-decay chains do not end with known nuclei 
making the method of a — a correlations insufficient for the unambiguous identification 
of isotopes. That these chains would end in unknown regions is easily understood when 
recognizing that with each successive a (or /5-decay), daughter nuclei move away from the 
regions of stability and spherical closed shells towards instability. It is reasonable to expect 
spontaneous fission (SF) to become the preferred decay mode over others. SF could then 
be considered to be a reliable signature for the formation of extremely heavy nuclei [ 33J. 
In such cases, and ideally suited to single atom production, chemical methods provide for 
excellent assignment criteria if the half-lives are within the range of detection i.e. around 
a few seconds, with cross-sections less than 100 pb following genetically established decays 
[ I34j . Generally, cold fusion experiments produce nuclei with shorter half-lives than those 
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synthesized by hot fusion. Hot fusion reactions are thus more suitable for chemical studies. 
Analogies may be drawn between the chemical properties of compounds of unknown Z 
with those of the chemical type of known elements allowing for specific assignments [ I35| . 
By adopting quick and sensitive techniques for separation, it has been possible to study 
Rf, Db, Sg and Bh (see [EE] for a complete review of chemistry). Furthermore, as a-decay 
may be considered to be the signature mechanism for the unique identification of each 
nuclide, genetic correlations are unambiguous. The detection of an a or a SF fragment 
is the only means for detecting a single atom following chemical separation which can 
be done with high efficiency. The method was used very effectively in confirming the 
a-decay grand-daughters for 277 112 synthesized at GSI, when 7 decay chains from Hs 
isotopes were unambiguously identified after chemical separation [EH]. First attempts to 
study the chemistry of Z = 112 were made by fusing natural U with a 48 Ca beam to form 
283 112 [ EE1 EH1 HOI • Due to various problems with the experimental set-up these results 
are still tentative and further experiments are underway [I41j. 

Where even shorter decay half-lives occur and as production cross - sections get smaller, 
still other techniques have been devised. Investigations wherein masses from evaporation 
residues (EVR's) are identified by exploiting the dependence of production cross-sections 
on excitation energy (thus defining the neutron yield) together with cross-bombardments 
in which the mass number of either the projectile or target is varied (thus changing the 
relative yield of resulting xn-evaporation channels) have been successfully carried out in 
the past. The methods find their use in both the identification of unknown nuclei as well 
as in the study of those with short SF decay half-lives [ 35]. This is best illustrated in 
the case of Z = 114 and Z = 116, where when re- visiting results from earlier experiments 
using more stringent methods, the Russian collaborations revised their previous conclu- 
sions resulting in re-assignments [ EU] for 286,287,288,289 114? 290,291,292,293 n6 [ - These 

experiments have been carried out using the Dubna Gas Filled Recoil Separator in 'hot 
fusion' reactions involving 233 > 238 U, 242 Pu and 248 Cm targets, collectively described in [ I42j . 
Other isotopes identified and discussed therein include 267 Rf, 271 Sg, 275 Hs, 279 ' 281 Ds, and 
282-285-Q2. j s expected that existing contradictions in the case of 287 114 and 283 112 will 
be resolved via continuing hot fusion studies along similar lines in the immediate future. 
Even these methods have their limitations where the use of cross-reactions and different 
decay channels become increasingly difficult with rapidly diminishing cross-sections for 
higher Z. 

It is clear that the direct determination of Z requires a stringent mass resolution a few 
mass units ~ 1 % to measure the difference between a xn and an axn evaporation channel 
for instance. The mass resolution should in fact be better than 1 amu at the level of 
300 amu (Am/m < 0.3 %) [ESI IB]- R is interesting that for all atoms heavier than U, 
masses have not been measured. This is due to the inherently small separation efficiency 
of elements in the range 89 < Z < 103. The identification of nuclei produced using in-flight 
recoil separators to date is based on the kinematic characteristics of the recoil products 
- the energy of the recoils and the emission angle with respect to beam direction. The 
necessity to retain the kinematic properties of the heavy elements being studied imposes 
conditions on target thickness, background conditions and relatively poor (factor of ~ 
6) suppression of reaction by-products [EBJ. New generation "isotope separator on-line" 
(ISOL) based machines such as 'MASHA' (Mass Analyzer for Superheavy Atoms) coming 
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up at JINR-Dubna will be able to directly determine the mass of separated atoms in the 
range 112 < Z < 120. The limitation on measurement is set by the shortest measurable 
half-life, T x / 2 ~ 0.5 s [ 46,. MASHA is expected to be operational shortly and this will 
open up new avenues for the study of the chemical properties of the SHE. 

3. Status of Theory 

As experimental methods have advanced considerably over the last couple of decades, 
several theoretical formulations have been developed and employed to describe the ground 
state properties of each of the nuclei along the a - decay chains including binding energies, 
Q a values etc. In spherical even-even nuclei, alpha transitions are expected to go from 
ground state to ground state and hence alpha decay energies (Q- values) serve as a sensitive 
signature for local changes in stability. This is not quite so clear in the case of deformed 
doubly even nuclei where transitions may originate from low lying rotational levels [ 14*7] 
or odd nuclei where transitions may be greatly hindered and spin isomerism may occur. 

Whereas observed Q- values may, in general, be measured quite accurately, uncertainties 
do exist in experimental half-lives. Cross-section limits were already less than 12 pb for the 
creation of Ds and Z=lll, for instance [EHHU3 and dropped to about 1 pb for the Element 
112 [H*5], resulting in poor statistics offered by single-atom count-rates. Interestingly, it is 
seen that a-decay half-lives are not as precise as one may think in microscopic calculations 
either: small variations in Q-values (to within a couple of hundred keV) dramatically 
affect theoretical estimates of half-lives. Given this imprecision in observed half-lives, 
they cannot be used in themselves for experimental assignments. However, a relatively 
long half-life in the presence of a, high Q a - value, serves as a strong indicator for elements Z 
> 100 [[515]. Furthermore, longer T a 's for even-even nuclei would imply higher hindrance 
factors, thereby serving to exclude the expected ground-to-ground transitions in such 
cases. Since the experimental presence of such hindrance factors could be considered as a 
property for assignment (or its exclusion) [E"5], theoretical predictions and comparisons 
will serve as a valuable guide. Finally, the qualitative trend of Q-values over say, neutron 
numbers, could serve as a good indication of the transition from deformation to sphericity 
as one traverses the neighborhood of N=162 shell and goes towards the spherical Z=114 
region. It is important that such arguments be kept in mind when comparing calculations 
with experiment. 

Several theoretical investigations have been carried out using the microscopic - macro- 
scopic (MicMac) method and the self consistent mean field in both the relativistic and 
non relativistic formalisms. The primary aim in early studies has been to predict the com- 
bination of neutron number (N) and proton number (Z) where spherical shell closure may 
occur. An "island of stability" had been predicted around the hypothetical doubly magic 
298 114 (N=184) about 30 years ago. More recently, nuclei in this vicinity are expected to 
be spherical or almost so with longer half-lives. Most theories do predict N=184 as being 
magic, however, there is no consensus on the location of the proton magic number due 
to differences in the treatment of the large Coulomb term and the spin-orbit interaction. 
MicMac models, which assume a prior knowledge about the densities and single particle 
potentials, include the Finite Range Droplet Model with folded Yukawa single particle po- 
tentials (FRDM+FY) [ 50. and the Yukawa plus Exponential model with Woods-Saxon 
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single particle potentials (YPE+WS) [ I51| E2], both of which confirm the prediction of 
298 114 18 4 as being the next spherical doubly-magic nucleus. Non relativistic microscopic 
models such as the Skyrme-Hartree-Fock-Bogoliubov method [ [53] where the spin-orbit 
term has to be manually introduced, predict Z=120 may be as probable as Z=114, indi- 
cating that magic shells in this region are isotope dependent [EH [55] . Such techniques 
tend to overestimate the splitting of levels due to the spin-orbit coupling which may effect 
predictions for shell closures. With the large density of single particle states which in 
turn characterizes this mass region, the SHE serve as a sensitive probe for distinguish- 
ing between the various theories that attempt to predict shell structure, especially when 
these models describe stable nuclei with comparable accuracy. Also, it has been known 
for some time that deformation effects are important to the understanding of stability 
in this region [ [53] • Bohr and Mottelson have observed that deformation may enhance 
stability [155]. 

Relativistic Mean Field theories (RMF) which incorporate the spin-orbit term naturally, 
[ EZl EHl EH1 E01 EH E21 ESI EU generally do better than the non-relativistic MF models. 
RMF has been very successful in describing ground state (g.s.) properties for nuclei 
spanning the entire periodic table, thereby claiming a 'global fit' [ |5JJ EU E5J. This 
method has also been applied successfully [ EH1 EZl EHJ EHl EDI IZH E2J ESI EH ESI to 
describe the ground states of these super-transactinides where both pairing effects and 
deformation play an important role. However, it is known that in all self-consistent models 
the occurrence of a spherical proton (neutron) shell closure with a given Z (N) can change 
with varying neutron number N (Z). Using the RMF formalism, neutron number N=162 
has been predicted [|55] to exhibit shell closure at around Z = 108 - 110. This is consistent 
with the findings of [E5J. In addition, the calculated shell corrections (see Tab. 5, page 
222 of [E5J) P eak at N = 166 for Z = 112 indicating a stability for Z = 112 around N = 
166, in agreement with experimental observations of the superheavy nucleus 277 112. 

In this work we present comprehensive and systematic calculations for experimentally 
observed chains of SHE in the region 109 < Z < 118. Binding energies, deformations, 
Q-values, radii and densities are calculated for the ground states and compared with the 
corresponding experimental values where available. The calculations agree well with the 
experiment. Most of the earlier investigations have been devoted to the description of the 
g.s. properties of SHE. Here, we focus on calculating their half lives. 

The calculated RMF densities are used to derive the projectile-daughter interaction 
energy using the Double Folding (DF) model with the density dependent M3Y (DDM3Y) 
effective nucleon-nucleon interaction. This in turn is used to calculate the a-decay half 
life of the parent nucleus using the WKB approximation [ [73 [75] . 

3.1. Relativistic Mean Field Theory 

The Relativistic Mean Field (RMF) theory [[5JJE3] is now established to be one of the 
most successful approach for the accurate description of nuclear properties. It starts with 
a Lagrangian describing the Dirac spinor nucleons interacting via exchange of mesons and 
the photon. The mesons considered here are: the isoscalar - scalar a, isoscalar - vector 
uo and isovector - vector p mesons. The a (u) meson produces long range attraction 
(short range repulsion), whereas the p meson is necessary for the isospin dependence of 
the nuclear properties. The photon, as usual, produces the Coulomb interaction. The 
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Lagrangian consists of free baryon and meson terms and the interaction terms. The 
Euler - Lagrange variational principle yields the equations of motion. In the mean field 
approximation, replacing the fields by their expectation values, one gets a set of non-linear 
coupled equations: 

1. The Dirac equation with potential terms involving meson and electromagnetic (e.m.) 
fields describing the nucleon dynamics 

2. A set of Klein-Gordon type equations with sources involving nucleonic currents and 
densities, for mesons and the photon. 

This set of equations, known as RMF equations is to be solved self - consistently. The 
pairing correlations, essential for the description of open shell nuclei, are incorporated 
either by simple BCS prescription, or self consistently through the Bogoliubov transfor- 
mation. The latter leads to the Relativistic Hartree Bogoliubov (RHB) equations. The 
RHB equations [ 64, read: 



h D - A 
—A* 





(1) 



Here, A is the Lagrange multiplier, Ek is the quasi-particle energy and U k and V k are the 
four dimensional Dirac spinors, normalized as: 

J (Ulu k , + V ] k V k ) = 5 kk , ; (2) 
ho is the usual RMF Dirac Hamiltonian: 

h D = -ia-V + P{M + g a a) + gju + g p T 3 p° + e 1 ~ n A° (3) 

Here M is the nucleon mass, a is the scalar field and u°, p% and A° are the Lorentz time 
like components of the respective meson and e.m. fields. These fields are to be determined 
self-consistently through the Klein-Gordon (KG) equations: 

{-V 2 + ml} a = - 9(7 p s - g 2 a 2 - g 3 a 3 (4) 

{-V 2 + ml} oo° = g wPv + 94 uj o3 (5) 

{-V 2 + m 2 p }pl = g p p 3 (6) 

-V 2 A° = ep c (7) 

Here, m a (g a ), m w (g u ) and m p (g p ) are the masses (coupling constants) of cr, uj and p 
fields respectively; g 2 is the coupling constant for the cubic self interaction terms for the 
a field; g 3 (g 4 ) is the coupling constant for the quartic self interaction term for the a {uj) 
field and e is the electronic charge. The sources (nuclear currents and densities) appearing 
in the above Klein-Gordon equations involve super spinors (U {V)): 

Ps = E V k h°V k , (8) 
E k > o 

pv = E ( 9 ) 

E k > 
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p 3 = E vhVk, (io) 

s fc > 

pc = E vtfi^v* (ii) 

s fc > 

In practice, the sum in Eqs. (jHl - ITT]) is taken over the positive energy states (no-sea 
approximation) . 

The RHB equations have two distinct parts: the self consistent field (/to) that describes 
the long range particle-hole correlations and the pairing field (A) that accounts for the 
correlations in the particle-particle (pp) channel. The pairing field A is expressed in 
terms of the matrix elements of the two body nuclear potential V pp in the pp-channel and 
the pairing tensor involving the super-spinors (U,V). In the case of the constant gap, 
A (= A) becomes diagonal resulting in the BCS type expressions for the occupation 
probabilities (v 2 ) [EH IE]: 



2 1 



1 efc - A 



\e k - A) + A2 



'121 



€k being the energy of the single particle state k and the Lagrange multiplier A (fermi 
energy) is to be determined through the BCS number equation. As a result, the RHB 
equations (Eq. |TJ) reduce to the RMF equations with a constant gap. 

A reliable and satisfactory derivation of V pp is not yet achieved in RMF [ |MJ [77] . 
Therefore, in practice, it is customary to adopt a phenomenological approach while solving 
the RHB equations. Therefore, one often uses for V pp , the finite range Gogny-DlS 
interaction, which is known to have the right pairing content and is given by: 



V(n,r 2 ) = E e ~ {(ri ~ r2)/M2 (W^ + BiP' 7 - HiP T - MiP°P T ). (13) 

=1,2 



where, /ij, Wi, B>i, Hi & Mj (i—1,2) are parameters of the interaction. In the case of the 
constant gap approximation, the required gap parameters are fixed so as to reproduce the 
corresponding Gogny D1S pairing energies. 

Next we present and discuss the calculated ground state properties of the nuclei be- 
longing to a - decay chains of superheavy nuclei. 



4. Results and Discussions 

4.1. Ground State Properties 

The RMF / RHB equations are solved either using the basis expansion technique, or 
in coordinate space. Explicit numerical calculations require the following inputs: 

1) The parameters appearing in the Lagrangian and 

2) pairing gaps (or V pp ) 

The output includes Dirac spinors, single particle (quasi-particle) energies, occupancies 
and the mesonic and e.m. fields. These in turn are used to obtain total binding energies, 
radii, densities, deformations etc. 

A number of Lagrangian parameter sets exist in the literature. Here, we employ the 
most frequently used parameter set NL3 [ ISO] . For comparison, we also use the NL-SV1 
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parameter set [ IHT] in some cases. These sets are listed in Tab. (1). In NL-SV1, the 
coupling constant (#4) for the quartic self - interaction term of the isoscalar - vector uo 
field is also included, and its value is 41.010. This inclusion improves the equation of 
state of nuclear matter. It is seen that results obtained by using the other Lagrangian 
parameter sets (e.g. NL1 [EHJUn], NL-SH exhibit similar systematics. Therefore, 

the inferences made and the conclusions drawn here are expected to hold for the other 
parameter sets as well. 

In this work we solve the RMF equations in the axially symmetric deformed oscillator 
basis. The resulting quantities (binding energies, deformation parameters, matter radii, 
etc.) are denoted by NL3 / NL-SV1. The required (spherical) pairing gaps are tuned 
to reproduce the pairing energy obtained from RHB calculations using the Gogny D1S 
interaction in the pairing channel. 

Table 1 

The different sets of Lagrangian parameters commonly used in the RMF / RHB calcu- 
lations. The masses are in MeV. All the coupling constants are dimensionless, except g<i 
which is expressed in terms of frrT 1 . 





M 


m a 




m p 


9a 


9u 


9 P 


92 


93 


NL3 


939 


508.194 


782.501 


763 


10.217 


12.868 


4.474 


-10.431 


-28.885 


NL-SV1 


939 


510.035 


783 


763 


10.125 


12.727 


4.492 


-9.241 


-15.388 



We now present and discuss the calculated ground state properties. The numerical 
values of the calculated ground state properties and the decay half lives are presented in 
Tab. H 

4.1.1. Binding Energies 

The calculated and the corresponding extrapolated (Audi2003) [ ISSj binding energies 
for the nuclei appearing in some of the a - decay sequences of observed superheavy nuclei 
are presented in Fig. (JTJ). The DEF results obtained using NL3 (NL-SV1) parameter set 
are denoted by NL3 (NL-SV1). It can be seen that both the calculations are in agreement 
with each other, and also with the corresponding Audi2003 values. However, it should be 
noted that the NL3 parameter set always yields a slightly deeper solution in comparison 
with NL-SV1. Further, both these calculations are found to yield consistently larger 
binding energies in comparison with the corresponding Audi2003 values. The differences, 
however, are small (about 5 parts in 2000). 

4.1.2. Deformations 

The quadrupole deformation parameters {(3) are obtained from the calculated point 
neutron (Q n ) and proton (Q p ) quadrupole moments through: 



Q = Q n +Q P = J^^ARlp 
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with R Q = 1.2A 1 / 3 (fm). The calculated (3 for the nuclei considered are shown in Fig. 
H2D, along with the corresponding Moller - Nix (MN) values [ 84j. The NL3 and NL-SV1 
parameter sets yield similar deformations. These graphs reveal a number of additional 
interesting features. First, most of the nuclei investigated here turn out to be prolate in 
shape. Very few of these are spherical. These findings agree with those reported by Moller 
and Nix [ 84J with some exceptions. Further, highly deformed solutions do appear for some 
very heavy nuclei with both NL3 and NL-SV1 parameter sets. We note that apart from 
these highly deformed solutions (which happen to be the lowest), nearby solutions also 
exist with smaller values of (3. These solutions are indicated by open / filled boxes in the 
respective figures for the deformation parameter (3. It has been shown in the literature [ 
E5] that with the incorporation of additional constraints, e.g. the octupole, these highly 
deformed solutions disappear. We shall ignore these solutions, and not discuss them here. 
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Figure 1. The calculated and the Audi2003 [IHH] binding energies for the nuclei belonging 
to a - decay chains of different superheavy nuclei. 
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Figure 2. The quadrupole deformation parameters. The corresponding Moller - Nix values 
[ ISlj are also indicated for comparison. Highly deformed solutions are indicated by open 
(filled) squares for NL-SVl (NL3) results. 
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Figure 3. Neutron pairing energies for even-even nuclei in the region 108<Z<118. 




Figure 4. Single particle states for Z = 116 isotopes. 



12 



4.1.3. Pairing Energy 

One of the characteristic features of a closed shell is the existence of a sharp maximum 
in the pairing energy curve. Here, we investigate the calculated D1S pairing energies for 
even-even nuclei in the isotopic chains of Z = 108 to Z = 118. For this purpose, the 
calculated neutron pairing energies corresponding to spherical solutions are shown in Fig. 
(J3J). This figure reveals that the pairing energy for Z = 108 at neutron number (N) 184 
is zero, indicating that this neutron number may indicate shell closure for Z = 108. As 
Z increases, the pairing energy for this neutron number also increases, while the position 
of the peak is unchanged. Beyond Z = 116, however, the peak disappears, which may 
indicate loss of 'magicity'. A similar 'peak' is observed at neutron numbers 172 and 164, 
where, the peak exists for all the elements considered. Further, for Z = 116 and 118, the 
pairing energy turns out to be zero. For elements lighter than Z=116, the pairing energy 
is deeper (i.e. negative) and turns out to be deepest at Z = 108. This analysis, therefore, 
indicates that the existence of such 'magicity' depends on the combination of both proton 
and neutron numbers rather than on either one alone. This conclusion is consistent with 
the findings of [tBUj. 

4.1.4. Single Particle States 

Magicity can be associated with the existence of a large gap between the last occupied 
state and the first unoccupied state. To examine this in detail, we plot the calculated 
(NL3) single particle states for isotopes of element 116 in Fig. The fermi levels are 
joined by a continuous solid line. A gap is evident at neutron number 172 ( 288 116). This 
observation is consistent with the plot of pairing energies where, at neutron number 172, 
the pairing energy has a sharp maximum around MeV. Furthermore, the isotopes of 
Z = 116 turn out to be spherical or nearly so, in the neighborhood of neutron number 
172 ( 288 116). This fact lends support to the conclusion that for the Z=116 isotopes, shell 
closure may be indicated at neutron number 172. 

To convey the structure of the single particle levels near the Fermi surface, we display 
the calculated (NL3) levels in Figs. (J5HZJ) for the nuclei appearing in the observed a - 
decay chains of SHE. 



Table 2: The calculated ground state properties of SHE. The Q values and decay half 
lives are also indicated. Here, Ex+WKB represents the half lives obtained by using the 
experimental Q values in WKB. 







BE/ A (MeV) 


P 


r a 


Q a (MeV) 


lo gio T «( s ) 


A 


Z 


NL3 


Audi03 


NL3 


MN 


(fm) 


NL3 


Expt. 
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BE/A (MeV) 
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Z 


NL3 


Audi03 


NL3 


MN 


(fm) 


NL3 


Expt. 


NL3 


Ex+WKB 


Expt. 


260 


105 


-7.38 
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0.24 
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-1.19 
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-0.24 
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0.97 
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-7.43 


-7.40 


0.28 


0.24 
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1.82 
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1.82 
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-7.43 
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257 


102 
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1.57 
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-7.42 


0.29 


0.24 
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101 
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249 


100 


-7.49 


-7.46 


0.29 


0.23 
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4.1.5. Q Values 

The Q value of a parent nucleus against a - decay is simply the difference between the 
binding energy of the parent nucleus and the sum of the binding energies of both the 
resulting daughter and a particle. Here, we present Q-values for the NL3 and NL-SV1 
parameter sets in Fig. (|Sjl. Both calculations are found to be in agreement with each 
other, and also with the corresponding experimental values. The calculations deviate 
from experiment. At some places, the maximum departure is of the order of about 1 
MeV, except for the decay chain on Z = 115, where the deviation is about 2 MeV in a few 
cases. This level of agreement in indeed gratifying in view of the fact that the Q value is 
the difference between large numbers. A small error in even one of them, could affect the 
Q value substantially. 

4.1.6. Matter Radii 

The root mean squared (rms) matter radii (r m ) are obtained from the rms proton and 
neutron radii through: 

„ Zrl + Nr 2 n 

2 p n / 14 x 

Z + N { ' 

The calculated r m values are shown in Fig. Q. The r m values are found to be varying 
monotonically with mass number. The radius parameter r Q extracted using the conven- 
tional relation: 

r m = r A V * (15) 

is found to be a constant (r = 0.9639 ± 0.0005) for all the nuclei appearing in the a - 
decay chains of SHE, shown in figure 
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Figure 5. Single particle level structure near the Fermi level 
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Figure 6. Single particle level structure near the Fermi level 
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Figure 7. Single particle level structure near the Fermi level 
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Figure 8. The Q values against a - decay. The corresponding experimental values are 
also shown, where available. 
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4.2. Half Lives 

4.2.1. Double Folding Model and WKB analysis 

The basic idea behind the double folding model is to get a reasonable nucleus-nucleus 
potential, knowing some nucleon-nucleon interaction (e.g. M3Y interaction, etc.). In 
general, the double folding potential comprises of the direct and the exchange terms. The 
direct term contains the direct nucleon - nucleon matrix elements, whereas, the exchange 
term, as the name suggests, contains the exchange part. The latter is considerably more 
difficult to handle in practice. Thus, for some of the applications, the exchange term is 
simulated by a delta function pseudo-potential, with some density dependence. 

In the present work, the double folding prescription to obtain the nucleus - nucleus 
potential has been used (refer to Fig. (jllj) for the geometry of the problem). The M3Y 
effective nucleon - nucleon interaction employed here is given by [ EEl EH EE] : 

—4s —2 5s 

v msy = 79m e 2134^— (16) 

The exchange effects are considered only through a delta function pseudo-potential [ IH7j : 
v pseudo = Joa{E)5(s) ; (17) 
where, the volume integral (Joo) is [ 87 : 

Joo = - 276(1 - °" 005E ) in MeV - fm 3 (18) 

In the present work, the energy dependence is ignored. Thus, the M3Y interaction with 
pseudo-potential becomes: 

—4s —2.5s 

v M3Y +P seudo = 79gg e 2 134— 276<5(s) in MeV (19) 

4s 2.5s 

The density dependence is supposed to compensate to some extent the higher order ex- 
change effects and the effects of the Pauli Blocking. Following the earlier work [EZI, it is 
assumed to be of the form: 

= C(l - P{E) P T){1 - P(E)pT) (20) 

where, C is the overall normalization constant, and is taken to be 1.0 in the present work. 
/3(E) is the energy dependent part of the density dependent term and is assumed to have 
a constant value, 1.6 [ 187] . 

With the density dependence and the M3Y with pseudo-potential, the assumed nucleon- 
nucleon interaction is now given by: 

v(s) = C (j999^ - 2134^ - 2765(s)) (1 - (3{E)pT){^ ~ (3{E)pT) (21) 

In these expressions, (3(E) is 1.6 and s is equal to f pt (refer to the figure). 
The total double folding potential between the nucleus - nucleus system is: 

Vpt(R) = J PP (f p )p T (f t )v(f p -f t + R)dr p dr t (22) 
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Figure 11. Geometry of the Double Folding prescription. 



where, v is as defined above (Eq. (|2ip). In the actual calculations, we evaluate the six 
dimensional integral in the above equation by transforming it into the momentum space. 
The details are easy to work out and are straightforward. 

In order to evaluate the half life time of the nucleus against a - decay, one needs to 
know the nuclear potential, the Coulomb potential and the energy of the alpha particle, 
which in turn is obtained by the Q value and the zero point energy of the oscillation of 
the alpha particle in the potential well. The point to be noted here is, in the present 
case, it is assumed that the alpha particle is already formed in the parent nucleus, so 
that, its motion can be simulated by assuming that it is moving in an average potential 
well formed by the daughter - a system. The Coulomb potential required in this case, is 
obtained by folding the point proton densities of target and projectile with the Coulomb 
interaction. 

Here the decay half life is calculated using the WKB approximation. The half life time 
of the nucleus (parent) against the a - decay is then given by: 

T 1/2 = ^ (l + e K ) (23) 

where, within the WKB approximation, the action integral is given by: 
2 r R b 

IRa 



K = -J R b {2n(E(R)-Q)} 1/2 dR, (24) 
and 



E(R)=V PT (R) + Vc(R) + n 2 



^A 2 
2/i r 2 

The third term represents the centrifugal barrier with A = I + 1/2, I being the orbital 
angular momentum. The R a and Rb are the lower and upper turning points respectively. 
These are determined through the requirement: 



E{R a ) = Q + E U = E{R b ) 
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In Eq. (J2"3|) . z/ is the assault frequency, given by: 

where, R is the radius of the parent, given by R = 1.2A 1 / 3 ; E is the energy of the alpha 
particle, corrected for recoil; M is mass of alpha particle, expressed in MeV. 
We now present and discuss the calculated half live values. 

4.2.2. Half Lives 
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Figure 12. Half lives against a - decay of superheavy nuclei. 



The calculated and the corresponding experimental half lives for the superheavy nuclei 
against a - decay are presented in Fig. (JT2J) - Here, we have not presented the results 



23 



for NL-SV1, since, we have already shown that the difference between NL3 and NL-SV1 
results is negligible. 

Results of two independent calculations have been presented in Fig. (J 12)) . In the first 
case, the calculated Q values are used in the WKB procedure. These results are denoted 
by NL3. In the second case, the experimental Q values are used in the WKB procedure, 
with the same microscopic potential, as in the first case. These results are denoted by 
Qexpt. + WKB in Fig. (|12j) . It is clearly seen that the use of experimental Q values in 
WKB approach reproduces the experimental values rather well, indicating the reliability 
of the microscopic nucleus - nucleus potential, and also validating the procedure that we 
follow in this work. The half lives obtained by using the calculated Q values, though 
similar to the experimental trend, differs from it quantitatively at places. This reflects 
hyper - sensitivity of the half lives on Q values. 

5. Summary and conclusions 

Extensive and systematic microscopic self consistent RMF calculations for the nuclei 
appearing in the a - decay chains of the known superheavy elements are presented. The 
RMF results with only a few (seven) fixed Lagrangian parameters are found to reproduce 
the corresponding experimental ground state properties remarkably well. This indeed, 
is gratifying. The predicted neutron shell closures N ~ 164, 172 and 184 are consistent 
with the recent experimental observations. The calculated root mean square radii closely 
follow the A 1 ^ law. 

The interaction potential between the a and daughter nucleus, obtained in the tpp 
approximation using the RMF densities and the nucleon - nucleon interaction M3Y is 
employed to estimate the a - decay half lives. The calculations qualitatively agree with 
the experiment. The use of experimental Q values brings the calculations closer to the 
experiment. Therefore, the calculated interaction potential is reliable and can be used 
with confidence in other reaction studies. 
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